Global Analysis of Fragmentation Functions for Eta Mesons 
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I. INTRODUCTION 



Fragmentation functions (FFs) are a key ingredient 
in the perturbative QCD (pQCD) description of pro- 
cesses with an observed hadron in the final-state. Sim- 
ilar to parton distribution functions (PDFs), which ac- 
count for the universal partonic structure of the inter- 
acting hadrons, FFs encode the non-perturbative details 
of the hadronization process When combined with 
the perturbatively calculable hard scattering cross sec- 
tions, FFs extend the ideas of factorization to a much 
wider class of processes ranging from hadron production 
in electron-positron annihilation to semi-inclusive deep- 
inelastic scattering (SIDIS) and hadron-hadron collisions 

1- 

Over the last years, our knowledge on FFs has im- 
proved dramatically [|| from first rough models of quark 
and gluon hadronization probabilities [J] to rather precise 
global analyses at next-to-leading order (NLO) accuracy 
of QCD, including estimates of uncertainties While 



'Electronic address: caidala@bnl.gov 



t Electronic address: ellingha@uni-mainz.de 
t Electronic address: sassot@df.uba.ar 
§ Electronic address: seelej@mit.edu] 

^Electronic address: marco@ribf.riken.jp address after Oct. 1 st : 
Physics Department, Brookhaven National Laboratory, Upton, 
NY 11973, USA 



the most accurate and clean information used to deter- 
mine FFs comes from single-inclusive electron-positron 
annihilation (SIA) into hadrons, such data do not al- 
low disentanglement of quark from anti-quark fragmen- 
tation and constrain the gluon fragmentation only weakly 
through scaling violations and sub-leading NLO correc- 
tions. Modern global QCD analyses [!, Hf utilize to the 
extent possible complementary measurements of hadron 
spectra obtained in SIDIS and hadron-hadron collisions 
to circumvent these shortcomings and to constrain FFs 
for all parton flavors individually. 

Besides the remarkable success of the pQCD approach 
in describing all the available data simultaneously, the 
picture emerging from such comprehensive studies reveals 
interesting and sometimes unexpected patterns between 
the FFs for different final-state hadrons. For instance, 
the strangeness-to-kaon fragmentation function obtained 
in Ref. [5| is considerably larger than those assumed pre- 
viously in analyses of SIA data alone Q. This has a 
considerable impact on the extraction of the amount of 
strangeness polarization in the nucleon [Toj from SIDIS 
data, which in turn is linked to the fundamental question 
of how the spin of the nucleon is composed of intrinsic 
spins and orbital angular momenta of quarks and gluons. 

Current analyses of FFs comprise pions, kaons, pro- 
tons [5-8], and lambdas 0, EH as final-state hadrons. 
In this respect, FFs are a much more versatile tool to 
explore non-perturbative aspects of QCD than PDFs 
where studies are mainly restricted to protons [HI, [l3| . 
In the following, we extend the global QCD analyses of 



2 



FFs at NLO accuracy as described in Refs. to eta 

mesons and estimate the respective uncertainties with 
the Lagrange multiplier method 0, [HI [3]. We ob- 
tain a parametrization from experimental data for single- 
inclusive eta meson production in SIA at various center- 
of-mass system (c.m.s.) energies y/S and proton-proton 
collisions at BNL-RHIC in a wide range of transverse 
momenta px- We note two earlier determinations of 
eta FFs in Refs. and [l6| which are based on nor- 
malizations taken from a Monte Carlo event generator 
and SU(3) model estimates, respectively. In both cases, 
parametrizations are not available. 

The newly obtained FFs provide fresh insight into the 
hadronization process by comparing to FFs for other 
hadrons. In particular, the peculiar wave function of the 
eta, \rj) ~ \uu + dd — 2ss), with all light quarks and 
anti-quarks being present, may reveal new patterns be- 
tween FFs for different partons and hadrons. The sim- 
ilar mass range of kaons and etas, m^o ~ 497.6 MeV 
and m n ~ 547.9 MeV, respectively, and the presence of 
strange quarks in both wave functions makes comparisons 
between the FFs for these mesons especially relevant. Of 
specific interest is also the apparently universal ratio of 
eta to neutral pion yields for pr > 2 GeV in hadron- 
hadrons collisions across a wide range of c.m.s. energies, 
see, e.g., Ref. [l7|, and how this is compatible with the 
extracted eta and pion FFs. 

In addition, the availability of eta FFs permits for the 
first time NLO pQCD calculations of double-spin asym- 
metries for single-inclusive eta meson production at high 
Pt which have been measured at RHIC [l8| recently. 
Such calculations are of topical interest for global QCD 
analyses of the spin structure of the nucleon |10j . Finally, 
the set of eta FFs also provides the baseline for studies of 
possible modifications in a nuclear medium fl9l 20| . for 
instance, in deuteron-gold collisions at RHIC [17|. 

The remainder of the paper is organized as follows: 
next, we give a brief outline of the analysis. In Sec. IIIII 
we present the results for the eta FFs, compare to data, 
and discuss our estimates of uncertainties. We conclude 
in Sec. HVl 



II. OUTLINE OF THE ANALYSIS 

A. Technical framework and parametrization 

The pQCD framework at NLO accuracy for the scale 
evolution of FFs plj and single-inclusive hadron produc- 
tion cross sections in SIA [22[ and hadron-hadron colli- 
sions [23| has been in place for quite some time and does 
not need to be repeated here. Likewise, the global QCD 
analysis of the eta FFs itself follows closely the methods 
outlined in a corresponding fit of pion and kaon FFs in 
Ref. [H], where all the details can be found. As in Q 
we use the Mellin technique as described in [To|, l24j to 
implement all NLO expressions. Here, we highlight the 
differences to similar analyses of pion and kaon FFs and 



discuss their consequences for our choice of the functional 
form parameterizing the FFs of the eta meson. 

As compared to lighter hadrons, in particular pions, 
data with identified eta mesons are less abundant and less 
precise. Most noticeable is the lack of any experimental 
information from SIDIS so far, which provided the most 
important constraints on the separation of contributions 
from Ujd, and s (anti-)quarks fragmenting into pions and 
kaons Q. Since no flavor-tagged data exist for SIA either, 
it is inevitable that a fit for eta FFs has considerably less 
discriminating power. Hence, instead of extracting the 
FFs for the light quarks and anti-quarks individually, we 
parametrize the flavor singlet combination at an input 
scale of no = 1 GeV, assuming that all FFs are equal, 
i.e., Dl = Dl = D[\ = D\ = D 7 ] = D|. We use the 
same flexible functional form as in Ref. Q with five fit 
parameters, 

DV(z, Ho) = N iZ a >(l- z) ft [1 + 7i (l - z) s >] x 

1 

B[2 + a,, A + 1] + ^B[2 + a u (3 t + 5 t + 1] ' (1) 

where z is the fraction of the four-momentum of the 
parton taken by the eta meson and i = u,u,d,d,s,s. 
B[a, b] denotes the Euler beta function with a and b cho- 
sen such that Ni is normalized to the second moment 
/p zDV{ z , ho) dz of the FFs. 

Although the assumption of equal light quark FFs 
seems to be rather restrictive at first, such an ansatz can 
be anticipated in view of the wave function of the eta me- 
son. One might expect a difference between strange and 
non-strange FFs though due to the larger mass of strange 
quarks, i.e., that the hadronization of u or d quarks is 
somewhat less likely as they need to pick up an ss pair 
from the vacuum to form the eta. Indeed, a "strangeness 
suppression" is found for kaon FFs [|| leading, for in- 
stance, to Dy > . In case of the eta wave function 
one can argue, however, that also a fragmenting s quark 
needs to pick up an ss pair from the vacuum. Never- 
theless, we have explicitly checked that the introduction 
of a second independent parameterization like in (JTJ to 
discriminate between the strange and non-strange FFs, 
does not improve the quality of the fit to the currently 
available data. Clearly, SIDIS data would be required to 
further refine our assumptions in the light quark sector 
in the future. 

The gluon-to-eta fragmentation is mainly con- 
strained by data from RHIC rather than scaling viola- 
tions in SIA. As for pion and kaon FFs in we find 
that a simplified functional form with 7 9 = in Eq. ([1]) 
provides enough flexibility to accommodate all data. 

Turning to the fragmentation of heavy charm and bot- 
tom quarks into eta mesons, we face the problem that 
none of the available data sets constraints their contri- 
butions significantly. Here, the lack of any flavor-tagged 
data from SIA hurts most as hadron-hadron cross sec- 
tions at RHIC energies do not receive any noticeable 
contributions from heavy quark fragmentation. Intro- 
ducing independent FFs for charm and bottom at their 
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respective mass thresholds improves the overall quality of 
the fit but their parameters are essentially unconstrained. 
For this reason, we checked that taking the shape of the 
much better constrained charm and bottom FFs for pi- 
ons, kaons, protons, and residual charged hadrons from 
[El, El, but allowing for different normalizations, leads to 
fits of comparable quality with only two additional free 
parameters. 

The best fit is obtained for the charm and bottom 
FFs from an analysis of residual charged hadrons p, i.e., 
hadrons other than pions, kaons, and protons, and hence 
we use 



D2(z,m c ) = D*(z,m e ) = N C D] 
D%(z,m b ) = Dl(z,m b )=N b D 



res 
b 



(z,m c ) 
(z,m b ) 



(2) 



N c and N b denote the normalizations for the charm and 
bottom fragmentation probabilities at their respective 
initial scales, to be constrained by the fit to data. The 
parameters specifying the D r c e b s can be found in Tab. Ill 
of Ref . @ . The FFs in Eq. §2$ are included discontinu- 
ously as massless partons in the scale evolution of the FFs 
above their MS thresholds /i = m c _ b with m c — 1.43 GeV 
and m b = 4.3 GeV denoting the mass of the charm and 
bottom quark, respectively. 

In total, the parameters introduced in Eqs. ((T|) and ([2]) 
to describe the FFs of quarks and gluons into eta mesons 
add up to 10. They are determined by a standard x 2 
minimization for N = 140 data points, where 



N 

£ 



SE] 



(3) 



Ej represents the experimentally measured value of a 
given observable, 5Ej its associated uncertainty, and Tj 
is the corresponding theoretical estimate calculated at 
NLO accuracy for a given set of parameters in Eqs. (fTJ) 
and For the experimental uncertainties 8Ei we take 
the statistical and systematic errors in quadrature for the 
time being. 



control which combinations of quark FFs are probed. 
Only the CERN-LEP data taken on the Z resonance re- 
ceive significant contributions from charm and bottom 
FFs. 

Given that the range of applicability for FFs is limited 
to medium-to- large values of the energy fraction z, as 
discussed, e.g., in Ref. @, data points with z < 0.1 are 
excluded from the fit. Whenever the data set is expressed 
in terms of the scaled three-momentum of the eta meson, 
i.e., x p = 2p Tj /y/S, we convert it to the usual scaling 

variable z = x p //3, where (3 — Pr)/E n = <Jl — m^/E^. In 

addition to the cut z > 0.1, we also impose that /? > 0.9 
in order to avoid kinematic regions where mass effects 
become increasingly relevant. The cut on j3 mainly affects 
the data at low z from BABAR [3^ |. 

In case of single-inclusive eta meson production in 
hadron-hadron collisions, we include data sets from 
PHENIX at V5 = 200 GeV at mid-rapidity [13, HI in 
our global analysis. The overall scale uncertainty of 9.7% 
in the PHENIX measurement is not included in 5Ej in 
Eq. ([3]). All data points have a transverse momentum 
Pt of at least 2 GeV. As we shall demonstrate below, 
these data provide an invaluable constraint on the quark 
and gluon-to-eta fragmentation probabilities. In general, 
hadron collision data probe FFs at fairly large momen- 
tum fractions z > 0.5, see, e.g., Fig. 6 in Ref. [Hj], 
complementing the information available from SIA. The 
large range of px values covered by the recent PHENIX 
data [23], 2 < p T < 20 GeV, also helps to constrain FFs 
through scaling violations. 

As in other analyses of FFs 0, Q we do not include 
eta meson production data from hadron-hadron collision 
experiments at much lower c.m.s. energies, like Fermilab- 
E706 39]. It is known that theoretical calculations at 
NLO accuracy do not reproduce such data very well with- 
out invoking resummations of threshold logarithms to all 
orders in pQCD 0. 



III. RESULTS 



Data sets included in the fit 



A total of 15 data sets is included in our analysis 



25] 
28] 



> 10 GeV: HRS 
= 29 GeV, JADE ^ 

35 GeV, and ALEPH 
[HI at \/S = Mz = 
Preliminary results from BABAR 36] at 



We use all SIA data with 
and MARK II [26| at 
and CELLO [29[ at VS = 34 - 
[Hil, L3 [H [H, and OPAL 
91.2 GeV 

VS= 10.54 GeV are also taken into account. 

The availability of e + e~ data in approximately three 
different energy regions of y/S ~ 10, 30, and 90 GeV helps 
to constrain the gluon fragmentation function from scal- 
ing violations. Also, the appropriate electroweak charges 
in the inclusive process e + e~ — > (7, Z) — > rjX vary with 
energy, see, e.g., App. A of Ref. [ll| for details, and hence 



In this Section we discuss in detail the results of our 
global analysis of FFs for eta mesons at NLO accuracy 
of QCD. First, we shall present the parameters of the 
optimum fits describing the at the input scale. Next, 
we compare our fits to the data used in the analysis and 
give x 2 values for each individual set of data. Finally, 
we estimate the uncertainties in the extraction of the 
using the Lagrange multiplier technique and discuss the 
role of the different data sets in constraining the FFs. 



A. Optimum fit to data 

In Tab. H] we list the set of parameters specifying the 
optimum fit of eta FFs at NLO accuracy in Eqs. (p} and 
@ at our input scale /j,q = 1 GeV for the light quark 
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flavors and the gluon. Charm and bottom FFs are in- 
cluded at their mass threshold /i = m c and /i = rrib, 
respectively [4lj |. 

The data sets included in our global analysis, as dis- 
cussed in Sec. Ill Bt and the individual x 2 values are pre- 
sented in Tab. HU We note that the quoted number of 

TABLE I: Parameters describing the NLO FFs for eta mesons, 
D^(z,fio), in Eqs. fTJ and @ at the input scale no = 1 GeV. 
Inputs for the charm and bottom FFs refer to fio = m c and 
fio — nib, respectively. 



Flavor i 


Ni 


on 


Pi 


7< 


St 


u, u, d, d, s, s 


0.038 


1.372 


1.487 


2000.0 


34.03 


9 


0.070 


10.00 


9.260 








c, c 


1.051 










b,b 


0.664 











points and x 2 values are based only on fitted data, i.e., 
z > 0.1 and > 0.9 in SIA. 

As can be seen, for most sets of data their partial con- 
tribution to the x 2 °f the fit is typically of the order of 
the number of data points or even smaller. The most 
notable exceptions are the HRS [H] and ALEPH '02 [32| 
data, where a relatively small number of points have a 
significant x 2 , which in turn leads to total % 2 per de- 
gree of freedom (d.o.f.) of about 1.6 for the fit. We have 
checked that these more problematic sets of data could 
be removed from the fit without reducing its constraining 
power or changing the obtained D 1 / significantly. The re- 
sulting, fairly large x 2 /d.o.f. due to a few isolated data 
points is a common characteristic of all extractions of 
FFs made so far 0-0] for other hadron species. 

The overall excellent agreement of our fit with experi- 
mental results for inclusive eta meson production in SIA 
and the tension with the HRS and ALEPH '02 data is 
also illustrated in Fig. [TJ It is worth pointing out that 
both ALEPH '00 [H| and BABAR [if data are well re- 
produced for all momentum fractions z in spite of being 
at opposite ends of the c.m.s. energy range covered by 
experiments. 

Our fit compares very well with all data on high-p^ eta 
meson production in proton-proton collisions from RHIC 
[H|33]. The latest set of PHENIX data [H significantly 
extends the range in px at much reduced uncertainties 
and provides stringent constraints on the FFs as we shall 
demonstrate below. The normalization and trend of the 
data is nicely reproduced over a wide kinematical range 
as can be inferred from Figs. [2HU In each case, the in- 
variant cross section for pp r\X at \/~S — 200 GeV 
is computed at NLO accuracy, averaged over the pseu- 
dorapidity range of PHENIX, \n\ <: 0.35, and using the 
NLO set of PDFs from CTEQ [ij] along with the cor- 
responding value of a s . Throughout our analysis we 
choose the transverse momentum of the produced eta as 
both the factorization and the renormalization scale, i.e., 
fi f = fi r = p T . 



Since the cross sections drop over several orders of 
magnitude in the given range of px, we show also the 
ratio (data-theory)/theory in the lower panels of Figs. [5^ 
131 to facilitate the comparison between data and our fit. 
One notices the trend of the theoretical estimates to over- 
shoot the data near the lowest values of transverse mo- 
menta, pt — 2 GeV which indicates that the factorized 
pQCD approach starts to fail. Compared to pion pro- 
duction at central pseudorapidities, see Fig. 6 in Ref. 
the breakdown of pQCD sets in at somewhat higher pt 
as is expected due to the larger mass of the eta meson. 

The shaded bands in Figs. EE] are obtained with the 
Lagrange multiplier method, see Sec. lIII Bl below. applied 
to each data point. They correspond to the maximum 
variation of the invariant cross section computed with 
alternative sets of eta FFs consistent with an increase of 
Ax 2 = 1 or Ax 2 = 2% in the total x 2 of the best global 
fit to all SIA and pp data. 

In addition to the experimental uncertainties prop- 
agated to the extracted D] 1 , a large theoretical ambi- 
guity is associated with the choice of the factorization 
and renormalization scales used in the calculation of the 
pp — > -qX cross sections. These errors are much more siz- 
able than experimental ones and very similar to those 
estimated for pp — > irX in Fig. 6 of Ref. [5J]. As in 
the DSS analysis for pion and kaon FFs [f| the choice 
fif = fi r = Pt and fif — /i r — S in pp collisions and SIA, 
respectively, leads to a nice global description of all data 
sets with a common universal set of eta FFs. 

Next, we shall present an overview of the obtained FFs 
D^(z, Q) for different parton flavors i and compare them 
to FFs for other hadrons. The upper row of panels in 



TABLE II: Data used in the global analysis of eta FFs, the 
individual % 2 values for each set, and the total \ 2 °f the fit. 



Experiment 


data points 


x 2 




fitted 




BABAR [36] 


18 


8.1 


HRS [25] 


13 


51.6 


MARK II [26] 


7 


3.8 


JADE '85 [2JJ 


1 


9.6 


JADE '90 [28] 


3 


1.2 


CELLO [29] 


4 


1.1 


ALEPH '92 [30J 


8 


2.0 


ALEPH '00 [31] 


18 


22.0 


ALEPH '02 [32] 


5 


61.6 


L3 '92 [33] 


3 


5.1 


L3 '94 [34] 


8 


10.5 


OPAL [35] 


9 


9.0 


PHENIX 2 7 [17] 


12 


4.1 


PHENIX 3tt [17] 


6 


2.9 


PHENIX '06 [37] 


25 


13.3 


TOTAL: 


140 


205.9 



5 



10 J r 



10' 



10 



10 



10 



10 



- 1 1 r - 



1 do" 




FIG. 1: Comparison of our NLO results with the data sets for inclusive eta meson production in SIA used in the fit; see Tab. 



Fig. [5] shows the dependence of the FFs on the energy 
fraction z taken by the eta meson at a scale Q equal to 
the mass of the Z boson, i.e., Q = Mz- Recall that at our 
input scale Q = fio = 1 GeV we assume that = D 7 ^ = 
D v d = Dj = D'J = Dg, which is preserved under scale 
evolution. At such a large scale Q = Mz the heavy quark 
FFs are of similar size, which is not too surprising as mass 
effects are negligible, i.e., m c ^ <C Mz- The gluon-to-eta 
fragmentation function is slightly smaller but rises 
towards smaller values of z. Overall both the shape and 
the hierarchy between the different FFs is similar 
to those found, for instance, for pions; see Fig. 18 in 
with the exception of the "unfavored" strangeness- 
to-pion fragmentation function which is suppressed. In 
order to make the comparison to FFs for other hadrons 
more explicit, we show in the lower three rows of Fig. [5] 
the ratios of the obtained D^(z, M z ) to the FFs for pions, 
kaons, protons from the DSS analysis 0, @j. 

The eta and pion production yields are known to be 
consistent with a constant ratio of about a half in a wide 
range of c.m.s. energies in hadronic collisions for px > 
2 GeV, but the ratio varies from approximately 0.2 at z ~ 
0.1 to about 0.5 for z > 0.4 in SIA [ijj. It is interesting to 
see how these findings are reflected in the ratios of the eta 
and neutral pion FFs for the individual parton flavors. 
We find that -D^+tz/ ^Z+a follows closely the trend of 
the SIA data as is expected since gluon fragmentation 
enters only at NLO in the cross section calculations. For 
strangeness the rate of eta to pion FFs increases towards 
larger z because of the absence of strange quarks in the 
pion wave functions. 



Inclusive hadron production at small-to-medium val- 
ues of pt is known to be dominated by gluon fragmen- 
tation at relatively large values of momentum fraction z 
0, HI] largely independent of the c.m.s. energy y/S. In 
the relevant range of z, 0.4 < z < 0.6, the ratio D^/Dg 
resembles the constant ratio of roughly 0.5 found in the 
eta-to-pion production yields. Both at larger and smaller 
values of z the Dg 1 is suppressed with respect to Dg . In 
general, one should keep in mind that FFs always ap- 
pear in complicated convolution integrals in theoretical 
cross section calculations [22|, |23[ which complicates any 
comparison of cross section and fragmentation function 
ratios for different hadrons. 

The comparison to the DSS kaon FFs |5J is shown in 
the panels in third row of Fig. [5] Most remarkable is the 
ratio of the gluon FFs, which is approximately constant, 
D^j/D^ ~ 2, over a wide region in z but drops below 
one for z > 0.6 At large z 



identical to D^ +a , 



DZ+u tends to be almost 



while D v s+§ 



resembles only at 



low z. The latter result might be understood due to the 
absence of strangeness suppression for Df +S , whereas a 
fragmenting s quark needs to pick up an s quark from 
the vacuum to form the eta meson. It should be noted, 
however, that kaon FFs have considerably larger uncer- 
tainties than pion FFs Q which makes the comparisons 
less conclusive. 

This is even more true for the proton FFs @. Nev- 
ertheless, it is interesting to compare our to those 
for protons which is done in the lower panels of Fig. [S] 
As for kaons, we observe a rather flat behavior of the 
ratio DV/DP, which drops below one at larger values of 
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THIS FIT 
Ax = 1 uncertainty 
Ax = 2% uncertainty 
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(data - theory) / theory 



f+tiii k t 



T 



2 4 6 8 10 12 

P T [GeV] 

FIG. 2: Upper panel: comparison of our NLO result for single- 
inclusive eta production in pp collisions at y/S = 200 GeV 
with PHENIX data where the eta is identified in the decay 
77 — > 27 [lTj. Lower panel: the ratio (data-theory) /theory. 
The shaded bands correspond to alternative fits consistent 



with an increase of A\ 2 
the best fit (see text). 



1 or Ax = 2% in the total \ of 



z. The corresponding rates for light quark FFs show the 
opposite trend and rise towards z — > 1. 

Regarding the relative sizes of the fragmentation prob- 
abilities for light quarks and gluons into the different 
hadron species, we find that eta FFs are suppressed 
w.r.t. pion FFs (except for strangeness), are roughly sim- 
ilar to those for kaons, and larger than the proton FFs. 
This can be qualitatively understood from the hierarchy 
of the respective hadron masses. For z > 0.6, the lack of 
decisive constraints from data prevents one from drawing 
any conclusions in this kinematic region. 

As we have already discussed in Sec. Ill Al due to 
the lack of any flavor tagged SIA data sensitive to the 
hadronization of charm and bottom quarks into eta 
mesons, we adopted the same functional form as for 
the fragmentation into residual charged hadrons @, i.e., 
hadrons other than pions, kaons, and protons. The fit 
favors a charm fragmentation almost identical to that 
for the residual hadrons (N c = 1.058) and a somewhat 
reduced distribution for bottom fragmentation (JVj = 
0.664). At variance to what is found for light quarks 
and gluons, after evolution, D' +5 and ^ differ signifi- 
cantly in size and shape from their counterparts for pions, 
kaons, and protons as can be also inferred from Fig. [5] 
Future data are clearly needed here for any meaningful 
comparison. 



[mb / GeV ] 

+ PHENIX data 
(t|->3ji) 




P T [GeV] 

FIG. 3: Same as in Fig.[2]but now for the data where the eta 
is identified in the decay 77 — > 3tt [l7j . 



10 
10 

10 ' 
10 ' 
10 ' 
10 ' 
10 

10 ' 

-] 

10 



,3 T| 

E^r [mb/GeV 2 ] 
dp' 

f PHENIX '06 data 



THIS FIT 

I Ax =1 uncertainty 
Ax _ =2% uncertainty 
-I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I- 



0.5 - (data - theory) / theory 


-0.5 



10 



15 20 

P T [GeV] 



FIG. 4: Same as in Fig. H but now for the latest PHENIX 
data H3]. 

B. Estimates of uncertainties 

Given the relatively small number of data points avail- 
able for the determination of the as compared to 
global fits of pion, kaon, and proton FFs [1, @ , we refrain 
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FIG. 5: Upper panels: individual FFs for eta mesons zD^(z,Q 2 ) at Q 2 = M 2 Z for i = u + u, s + s, g, c + c, and b + b. Lower 
three rows of panels: ratios of our eta FFs to the ones for pions, kaons, and protons as obtained in the DSS analysis [1,0]. 



from performing a full-fledged error analysis. However, 
in order to get some idea of the uncertainties of the 
associated with experimental errors, how they propagate 
into observables, and the role of the different data sets in 
constraining the , we perform a brief study based on 
Lagrange multipliers HEJG3. 

This method relates the range of variation of a physical 
observable O dependent on FFs to the variation in the 
X 2 function used to judge the goodness of the fit. To this 
end, one minimizes the function 

$(A,{a i })= X 2 ({ai}) + AO({a i }) 



(4) 



with respect to the set of parameters {a.i} describing the 
FFs in Eqs. (P) and |(2]) for fixed values of A. Each of 
the Lagrange multipliers A is related to an observable 
0({dj}), and the choice A = corresponds to the opti- 
mum global fit. From a series of fits for different values 
of A one can map out the x 2 profile for any observable 
0({di}) free of the assumptions made in the traditional 
Hessian approach (42j. 

As a first example and following the DSS analyses [U, 
H[, we discuss the range of variation of the truncated 
second moments of the eta FFs, 



£i (Zmin, Q) 



zDV( z ,Q) dz, 



(5) 



for 



0.2 and Q = 5GeV around the values ob- 



tained in the optimum fit to data, £^ . In a LO approx- 
imation, the second moments f Q zD^(z,Q)dz represent 
the energy fraction of the parent parton of flavor i taken 
by the eta meson at a scale Q. The truncated moments 
in Eq. ([5]) discard the low-z contributions, which are not 
constrained by data and, more importantly, where the 
framework of FFs does not apply. In general, FFs enter 
calculations of cross sections as convolutions over a wide 
range of z, and, consequently, the ^(z m in, Q) give a first, 
rough idea of how uncertainties in the FFs will propagate 
to observables. 

The solid lines in Fig. [5] show the £^(z m in,<2) defined 
in Eq. ([5]) for i = u + u, g, c + c, and b + b against 
the corresponding increase A^ 2 in the total x 2 of the 
fit. The two horizontal lines indicate a A\ 2 of one unit 
and an increase by 2% which amounts to about 4 units 
in x 2 , see Tab. |TTJ The latter A\ 2 should give a more 
faithful estimate of the relevant uncertainties in global 
QCD analyses || H, E3, E3 than an increase by one unit. 

As can be seen, the truncated moment associ- 
ated with light quark FFs Z^ + - = D v ^ = D v s+ -, is 
constrained within a range of variation of approximately 
-20% aroun d the value computed with the best fit, as- 
suming a conservative increase in \ 2 by 2%. The es- 
timated uncertainties are considerably larger than the 
corresponding ones found for pion and kaon FFs, which 
are typically of the order of ±3% and ±10% for the light 
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FIG. 6: Profiles of x 2 for the NLO eta fragmentation fit as a function of the truncated second moments ^(zmin = 0.2, Q = 5 GeV) 
for different flavors (solid lines). In each case, the moments are normalized to the value £o they take in the best fit to data. 
The other lines indicate the partial contributions to Ax 2 of the individual data sets used in the fit. 



quark flavors [5|j, respectively, but closer to the ±20% ob- 
served for proton and anti-proton FFs jf|. For the trun- 
cated moment £^ of gluons shown in the upper right panel 
of Fig. [(J the range of uncertainty is slightly smaller than 
one found for light quarks and amounts to about ±15%. 
The allowed variations are larger for charm and bottom 
FFs as can be inferred from the lower row of plots in 
Fig. 1 

Apart from larger experimental uncertainties and the 
much smaller amount of SIA data for identified eta 
mesons, the lack of any information from SIDIS is partic- 
ularly responsible for the large range of variations found 
for the light quarks in Fig. [5] We recall that the miss- 
ing SIDIS data for produced eta mesons also forced us to 
assume that all light quark FFs are the same in Eq. ([1}. 
The additional ambiguities due to this assumption are 
not reflected in the \ 2 profiles shown in Fig. [51 The FFs 
for charm and bottom quarks into eta mesons suffer most 
from the lack of flavor tagged data in SIA. 

To further illuminate the role of the different data sets 
in constraining the we give also the partial contribu- 
tions to Ax 2 of the individual data sets from pp collisions 
and the combined SIA data in all panels of Fig. [51 Sur- 



prisingly, the light quark FFs are constrained best by 
the PHENIX pp data from run '06 and not by SIA data. 
SIA data alone would prefer a smaller value for 
by about 10%, strongly correlated to larger moments for 
charm and bottom fragmentation, but the minimum in 
the x 2 profile is much less pronounced and very shallow, 
resulting in rather sizable uncertainties. 

This unexpected result is most likely due to the fact 
that the SIA data from LEP experiments constrain 
mainly the flavor singlet combination, i.e., the sum of 
all quark flavors, including charm and bottom. Since 
there are no flavor tagged data available from SIA for 
eta mesons, the separation into contributions from light 
and heavy quark FFs is largely unconstrained by SIA 
data. Only the fairly precise data from BABAR at 
y/S ~ 10 GeV provide some guidance as they constrain a 
different combination of the light u, d, and s quark FFs 
weighted by the respective electric charges. Altogether, 
this seems to have a negative impact on the constraining 
power of the SIA data. For not too large values of pt, 
data obtained in pp collisions are in turn mainly sensitive 
to D^j but in a limited range of z, 0.4 < z < 0.6, as men- 
tioned above. Through the scale evolution, which couples 
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FIG. 7: Relative fractions of quarks and gluons fragmenting 
into observed eta meson in pp collisions at VS = 200 GeV 
and PHENIX kinematics [13 ■ 

quark and gluon FFs, these data provide a constraint on 
Cu+u- I n addition, the latest PHENIX data extend to 
a region of px where quark fragmentation becomes im- 
portant as well. To illustrate this quantitatively, Fig. [7] 
shows the relative fractions of quarks and gluons frag- 
menting into the observed eta meson as a function of pt 
in pp collisions for PHENIX kinematics. As can be seen, 
quark-to-eta FFs become dominant for pr > 10 GeV. 

The x 2 profile for the truncated moment of the gluon, 
is the result of an interplay between the PHENIX 
run '06 pp data and the SIA data sets which constrain 
the moment & towards smaller and larger values, re- 
spectively. This highlights the complementarity of the 
pp and SIA data. SIA data have an impact on £^ mainly 
through the scale evolution in the energy range from LEP 
to BABAR. In addition, SIA data provide information in 
the entire range of z, whereas the pp data constrain only 
the large z part of the truncated moment Conse- 
quently, the corresponding \ 2 profile for z m - ln = 0.4 or 
0.5 would be much more dominated by jm data. In gen- 
eral, the other data sets from PHENIX [17J do not have 
a significant impact on any of the truncated moments 
shown in Fig. [6] due to their limited precision and cov- 
ered kinematic range. 

Compared to pion and kaon FFs [5j, all \ 2 profiles 
in Fig. [5] are significantly less parabolic, which prevents 
one from using the Hessian method [42] for estimating 
uncertainties. More importantly, the shapes of the x 2 
profiles reflect the very limited experimental information 
presently available to extract eta FFs for all flavors reli- 
ably. Another indication in that direction are the differ- 
ent preferred minima for the values of the by the SIA 
and pp data, although tolerable within the large uncer- 
tainties. Our fit is still partially driven by the set of as- 
sumptions on the functional form of and relations among 



different FFs, which we are forced to impose in order to 
keep the number of free fit parameters at level such that 
they can be actually determined by data. Future mea- 
surements of eta production in SIA, pp collisions, and, in 
particular, SIDIS are clearly needed to test the assump- 
tions made in our analysis and to further constrain the 

The large variations found for the individual FFs in 
Fig. [6] are strongly correlated, and, therefore, their im- 
pact on uncertainty estimates might be significantly re- 
duced for certain observables. If, in addition, the ob- 
servable of interest is only sensitive to a limited range of 
hadron momentum fractions z, than the corresponding 
X 2 profile may assume a more parabolic shape. 

In order to illustrate this for a specific example, we 
compute the % 2 profiles related to variations in the theo- 
retical estimates of the single-inclusive production of eta 
mesons in pp collisions at PHENIX kinematics 37]. The 
results are shown in Fig. [5] for four different values of pr 
along with the individual contributions to A% 2 from the 
SIA and pp data sets. As anticipated, we find a rather 
different picture as compared to Fig. O with variations 
only ranging from 5 to 10% depending on the px value 
and tolerating Ax 2 /x 2 — 2%- The corresponding uncer- 
tainty bands are also plotted in Fig. 0] above for both 
Ax 2 = 1 and Ax 2 /x 2 — 2% and have been obtained for 
the other pp data from PHENIX ^ shown in Figs. [21 
and [3] as well. 

The uncertainties for pp — > rjX are smallest for inter- 
mediate pr values, where the latest PHENIX measure- 
ment [37l | is most precise and the three data sets [I?], H3] 
have maximum overlap, and increase towards either end 
of the pt range of the run '06 data. In particular at inter- 
mediate pt values, the main constraint comes from the 
PHENIX run '06 data, whereas SIA data become increas- 
ingly relevant at low pt ■ The previous pp measurements 
from PHENIX [H are limited to p T < H GeV and have 
considerably larger uncertainties and, hence, less impact 
on the fit. 



IV. CONCLUSIONS 

A first global QCD analysis of eta fragmentation func- 
tions at NLO accuracy has been presented based on the 
world data from electron-positron annihilation experi- 
ments and latest results from proton-proton collisions. 
The obtained parameterizations [4l| reproduce all data 
sets very well over a wide kinematic range. 

Even though the constraints imposed on the eta meson 
fragmentation functions by presently available data are 
significantly weaker than those for pions or kaons, the 
availability of eta FFs extends the applicability of the 
pQCD framework to new observables of topical interest. 
Among them are the double-spin asymmetry for eta pro- 
duction in longitudinally polarized proton-proton colli- 
sions at RHIC, eta meson production at the LHC, pos- 
sible medium modifications in the hadronization in the 
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FIG. 8: As in Fig. [6] but now as a function of the cross section for pp — > r/ X at different values of pr- The variations of the 
cross sections are normalized to the value ao obtained in the optimum fit to data. 



presence of a heavy nucleus, and predictions for future 
semi-inclusive lepton-nucleon scattering experiments. 

The obtained FFs still depend on certain assumptions, 
like SU(3) symmetry for the light quarks, dictated by 
the lack of data constraining the flavor separation suffi- 
ciently well. Compared to FFs for other hadrons they 
show interesting patterns of similarities and differences 
which can be further tested with future data. 
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